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ABSTRACT 

We present the results of a stellar population analysis of 30 Lyman alpha emitting galaxies (LAEs) at z ~ 0.3, 
previously discovered with the Galaxy Evolution Explorer (GALEX). With a few exceptions, we can accurately 
fit model spectral energy distributions to these objects, representing the first time this has been done for a large 
sample of LAEs at z < 3, a gap of ~ 8 Gyr in the history of the Universe. From the 26/30 LAEs which we can 
fit, we find an age and stellar mass range of 200 Myr - 10 Gyr and 10 9 - 10 11 M , respectively. These objects 
thus appear to be significantly older and more massive than LAEs at high-redshift. We also find that these 
LAEs show a mild trend towards higher metallicity than those at high redshift, as well as a tighter range of 
dust attenuation and interstellar medium geometry. These results suggest that low-redshift LAEs have evolved 
significantly from those at high redshift. 

Subject headings: galaxies: fundamental parameters - galaxies: evolution 



1. INTRODUCTION 

Over the past decade thousands of Lyman alpha emitting 
galaxies (LAEs) have been discovered from 2.4 < z < 6.96 
(e.g., Cowie & Hu 1998; Rhoads et al. 2000; Rhoads et al. 
2001; Gawiser et al. 2006a; Kashikawa et al. 2006; Ouchi et 
al. 2008). It was predicted that strong Lya emission would be 
a signpost of primitive galaxies in formation, and thus these 
objects could represent some of the first galaxies in the Uni- 
verse (Partridge & Peebles 1967). However, with the new 
results coming from stellar population modeling, the identi- 
fication of LAEs as primitive galaxies has come into doubt. 
Numerous studies have begun to detect dust in LAEs, which 
is an indicator that a prior generation of stars has lived and 
died, thus making these objects not primitive (Chary et al. 
2005; Lai et al. 2007; Pirzkal et al. 2007; Nilsson et al. 
2007; Kuiper et al. 2009; Yuma et al. 2009; Finkelstein 
et al. 2008, 2009). While the presence of dust was origi- 
nally proposed as the reason why LAEs were not detected 20 
years ago, it is also possible that the dust can help enhance the 
value of the Lya EW by selectively attenuating the continuum 
(Neufeld et al. 1991; Haiman & Spaans 1999; Hansen & Oh 
2006). In our previous work (Finkelstein et al. 2008, 2009), 
we demonstrated that dust enhancement from a clumpy inter- 
stellar medium (ISM) can help explain the spectral energy dis- 
tributions (SEDs) of a sample of z ~ 4.5 LAEs, with clumpy 
dust appearing to exist in 10 out of a sample of 14 LAEs, even 
though 12/14 of these LAEs were young (< 15 Myr). We also 
found plausible evidence for dust in every object in our sam- 
ple, although their derived metallicities were rather low (with 
half of the sample having Z < 0.2 Z ). This suggests that 
these objects are not primordial, as a previous generation of 
stars was likely responsible for the creation of the dust. 

While we have estimates for these fundamental physical pa- 
rameters of distant LAEs, it becomes interesting to study them 
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at lower redshift where we can observe them in greater de- 
tail. However, large samples of LAEs have only been studied 
down to z ~ 3.1 (Gawiser et al. 2006a; McLinden et al. 2009; 
Nakamura et al. 2009), with a few recent samples currently 
under analysis at z ~ 2 (Nilsson et al. 2008; Reddy et al. 
2008). Given the loss of CCD quantum efficiency at ~ 3600 
A (and the atmospheric ultraviolet (UV) cut-off at ~ 3100 A), 
detection of LAEs at z < 2 is not possible from the ground. 
However, LAEs could be detected using the space-based UV 
telescope, Galaxy Evolution Explorer (GALEX), at z < 1 . De- 
harveng et al. (2008) recently published the discovery of 96 
LAEs at 0.2 < z < 0.35 via GALEX spectroscopy, using data 
from the publicly available GALEX GR2 data release. We 
present the results from a stellar population modeling analysis 
of 30 of these objects, with an emphasis on their differences 
from high-redshift LAEs. 

Where applicable, we assume cosmology with il m = 0.3, 
fl A = 0.7 and Ho = 70 km s" 1 Mpc" 1 (c.f. Spergel et al. 2007). 
All magnitudes in this paper are listed in AB magnitudes (Oke 
& Gunn 1983). 



2. Lya GALAXY SAMPLE 

Deharveng et al. (2008; hereafter D08) recently published 
a GALEX catalog of low-redshift Lya emitting galaxies, pro- 
viding their positions, spectroscopic redshifts, Lya line fluxes 
and EWs, and GALEX far-UV (FUV) and near-UV (NUV) 
broadband magnitudes. To compare these objects to stellar 
population models, we require optical broadband fluxes. We 
thus examined a subset of these LAEs in the Extended Groth 
Strip (EGS) and the Extended Chandra Deep Field - South 
(ECDF-S), both of which have extensive multi-wavelength 
public datasets. 

The All-Wavelength Extended Groth Strip International 
Survey (AEGIS; Davis et al. 2007) covers the EGS from radio 
to X-ray wavelengths. We used optical data from the Canada- 
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FIG. 1 . — Cutouts of the positions of the 30 LAE; MUSYC V-band for the CDFS, and CFHTLS i'-band for the EGS. Large solid-line circles are centered on the 
GALEX positions of these objects as published in Deharveng et al. (2008), with radii of 2.5", comparable to the FWHM of the GALEX beam. Small dashed-line 
circles denote all objects in the CFHTLS or MUSYC optical catalogs. Smaller solid-line circles are centered on the object selected by our matching code to be 
the optical counterpart. In the majority of cases, the circles all line up, meaning that the GALEXbeam is centered on the optical counterpart. 



France-Hawaii Telescope Legacy Survey 1 (CFHTLS), which 
covers one deg 2 of sky (with seeing < 1")- The 5 a limit- 
ing AB magnitudes of these data are: u' (27.2), g' (27.5), r' 
(27.2), i' (27.0) and z' (26.0). D08 publish the positions of 
39 Lya galaxies in the EGS, of which 30 fall in the CFHTLS 
area. We then use positional matching to locate the optical 
counterparts to these 30 objects. One object had no nearby 
optical counterpart, so we exclude it from further analysis, 
leaving us with 29 objects. Given that GALEX has a spatial 
resolution of 4" and 5.6" FWHM in the FUV and NUV chan- 
nels, respectively, GALEX positions can have high uncertain- 
ties relative to typical ground-based seeing (~ 5")- To account 
for this, we use a search radius of 7", and then visually in- 
spect the results to be sure that we are selecting the correct 
optical counterpart. Figure 1 shows these results, displaying 
the CFHTLS i'-band image overplotted with large solid-line 
circles denoting the GALEX positions (with a 5" diameter), 
dashed-line circles denoting all objects from the CFHTLS i'- 
band selected catalog (2" diameter), and smaller solid-line 
circles encompassing the i'-band source that we selected as 
the optical counterpart. In most cases, we see both large and 
small solid-line circles centered around an object. However, 

1 Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Tele- 
scope (CFHT) which is operated by the National Research Council (NRC) 
of Canada, the Institut National des Science de l'Univers of the Centre Na- 
tional de la Recherche Scientifique (CNRS) of France, and the University of 
Hawaii. This work is based in part on data products produced at the Cana- 
dian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope 
Legacy Survey, a collaborative project of NRC and CNRS. 



in two cases the GALEX position was between two or more 
sources, and in these latter cases we removed these objects 
from further consideration, leaving us with 27 EGS LAEs. It 
is worth noting that some of the excluded objects could pos- 
sibly be mergers, thus triggering star-formation and creating 
Lya photons. If this is the case, then the Lya could be coming 
from one or both objects, however a detailed analysis of this 
possibility will await higher resolution FUV or Lya imaging. 

We performed a similar procedure in the ~ 30' x 
30' ECDF-S, using data from the Multiwavelength Survey by 
Yale-Chile (MUSYC; Gawiser et al. 2006b), which cover the 
Extended ECDF-S to a 5 a depth of 26.0, 26.9, 26.4, 26.4, 
24.6 and 23.6 AB mag in the U, B, V, R, I and z' filters, re- 
spectively. Out of the 15 GALEX selected LAEs, we find that 
four are covered by the MUSYC data. Two of these have un- 
ambiguous optical counterparts, while a third has one bright 
object at the center of the GALEX beam, with another faint 
object near the edge of the beam. The fourth object has noth- 
ing at the center of the GALEX beam, with two objects at the 
edge. We thus kept the first three objects, while excluding the 
fourth one. Thus our final sample consists of 30 LAEs, 27 in 
the EGS and 3 in the ECDF-S. 

2.1. Stellar Population Modeling 

Our motivation for this study was to learn about the stel- 
lar population ages and masses of low-redshift LAEs. To do 
this, we use the stellar population modeling code of Bruzual 
& Chariot (2003; BC03), combined with the method from our 
previous studies (Finkelstein et al. 2007, 2008, 2009). Briefly, 
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TABLE 1 

Best-Fit Single Population Models of z ~ 0.3 LAEs. 



IAU Format Name 


Label 


tpoi 




Mass 




Z 


TSFH 


A 1200 


q 


EW 




z* 






(Myr) 


(10 9 M 


3) 


(Z ) 


(yr) 


(mag) 




(A) 






J03321 1.92-280130.0 


CDFS1 


6000.0 ± 


2695.5 


6.91 ± 


0.47 


0.005 


4xKF 


0.20 ± 0.37 


0.25 ± 1.38 


155 


36.11 


0.22 


J033256.66-275316.5 


CDFS2 


45.0 ± 


0.0 


4.56 ± 


0.15 


2.5 


10'' 


2.00 ± 0.00 


0.00 ± 0.13 


232 


29.73 


0.37 


J033307. 32-274432.6 


CDFS3 


3000.0 ± 


250.0 


107.70 ± 


9.46 


0.4 


10 s 


1.25 ± 0.35 


2.00 ± 0.50 


35 


12.52 


0.22 


J142041. 17+530650.2 


EGS1 


453.5 ± 


3273.3 


55.85 ± 


88.35 


1.0 


10' 


4.50 ± 2.07 


0.75 ± 1.38 


240 


16.76 


0.20 


J141743.43+522805.9 


EGS2 


4.0 ± 


0.0 


0.012 ± 


0.002 


0.2 


10- 


0.40 ± 0.05 


0.00 ± 1.50 


221 


371.17 


0.21 


J141814.1 1+522343.8 


EGS3 


404.2 ± 


719.9 


1.99 ± 


0.67 


0.4 


10' 


1.75 ± 0.64 


0.25 ± 0.75 


31 1 


9.50 


0.21 


J141710.64+-53 1 153.2 


EGS4 


5000.0 ± 


1000.0 


183.40 ± 131.80 


0.2 


10 s 


0.80 ± 1.10 


0.75 ± 0.25 


147 


16.34 


0.21 


J 142226.54+523027.0 


EGS5 


2500.0 ± 


2745.6 


^ Q1 -t- 


U.o 1 


0.02 


10 5 


0.20 ± 0.64 


0.25 ± 0.88 


135 


3.22 


0.22 


J 141 822.42+52 1823.6 


EGS6 


321.0 ± 


93.9 


1.11 ± 


0.20 


0.4 


10' 


0.70 ± 0.37 


0.00 ± 1.00 


179 


5.19 


0.24 


J141805. 12+524507.1 


EGS7 


286.1 ± 


15.6 


1.40 ± 


0.14 


0.4 


10' 


0.20 ± 0.23 


1.50 ± 0.88 


87 


7.58 


0.25 


J141848. 25+521756.1 


EGS8 


360.2 ± 


41.6 


1.28 ± 


0.22 


0.2 


10' 


0.40 ± 0.37 


0.25 ± 0.88 


136 


2.71 


0.24 


J141654.27+522440.0 


EGS9 


2500.0 ± 


4625.0 


24.95 ± 


18.19 


0.2 


10 s 


1.00 ± 0.62 


1.00 ± 0.63 


1 15 


2.98 


0.25 


J142043. 29+524307. 8 


EGS10 


227.3 ± 


206.6 


1.20 ± 


0.20 


0.4 


10' 


0.60 ± 0.32 


1.00 ± 1.38 


96 


1.08 


0.25 


J 1419 14.72+522326.4 


EGS11 


321.0 ± 


1178.3 


3.47 ± 


3.12 


0.2 


10' 


1.50 ± 1.01 


1.00 ± 2.13 


103 


1.52 


0.26 


J142044.64+525006.5 


EGS12 


286.1 ± 


856.9 


3.47 ± 


0.79 


0.2 


10' 


1.50 ± 0.62 


1.00 ± 0.75 


103 


2.63 


0.26 


J 1417 12. 18+523556.6 


EGS13 


1800.0 ± 


2106.9 


6.63 ± 


1.21 


0.02 


10' 


0.60 ± 0.67 


2.00 ± 1.00 


69 


3.46 


0.26 


J142124.54+523920.1 


EGS14 


3000.0 ± 


1000.0 


37.39 ± 


6.33 


0.2 


10" 


1.75 ± 1.06 


1.00 ± 0.63 


1 15 


1.57 


0.26 


J141855. 63+524927. 7 


EGS15 


4000.0 ± 


1500.0 


64.19 ± 


10.83 


0.4 


10" 


0.30 ± 0.60 


3.00 ± 1.75 


62 


5.40 


0.27 


J 141952.94+522100.2 


EGS16 


3000.0 ± 


4361.0 


62.34 ± 


7.71 


0.4 


10" 


1.50 ± 1.20 


1.50 ± 1.00 


53 


4.78 


0.27 


J141907. 21+522108. 3 


EGS17 


5000.0 ± 


2179.7 


151.77 ± 


31.44 


0.2 


10 s 


0.50 ± 0.97 


1.50 ± 0.63 


98 


5.66 


0.27 


J141930.24+530204.8 


EGS18 


286.1 ± 


615.7 


2.22 ± 


0.40 


0.2 


I0 7 


1.00 ±0.51 


1.50 ± 1.13 


66 


2.22 


0.27 


J141709.87+530511.3 


EGS19 


360.2 ± 


1339.5 


5.37 ± 


2.97 


0.02 


I0 7 


1.75 ±0.83 


0.50 ± 1.25 


225 


4.43 


0.27 


J141936.54+530936.2 


EGS20 


321.0 ± 


37.0 


1.73 ± 


0.27 


0.02 


I0 7 


0.60 ± 0.32 


1.00 ± 1.38 


103 


2.11 


0.27 


J141937.64+523024.8 


EGS21 


1800.0 ± 


4492.4 


14.72 ± 


2.54 


0.2 


10 s 


1.25 ±0.67 


1.50 ± 1.75 


66 


1.55 


0.29 


J141649.90+525019.9 


EGS22 


2500.0 ± 


1856.9 


10.12 ± 


3.52 


0.02 


10 s 


0.20 ± 0.83 


5.00 ±2.13 


58 


2.35 


0.29 


J141855.91+525931.8 


EGS23 


202.6 ± 


615.7 


1.35 ± 


0.32 


0.2 


l() 7 


1.00 ±0.55 


1.00 ±0.63 


104 


2.72 


0.29 


J142208. 10+525225.8 


EGS24 


508.8 ± 


27.7 


37.99 ± 


6.58 


0.2 


I0 7 


2.50 ± 0.46 


0.25 ±0.13 


566 


6.27 


0.31 


J141616.62+530133.3 


EGS25 


5.0 ± 


317.8 


0.36 ± 


1.14 


1.0 


10 5 


3.50 ± 1.20 


0.75 ± 0.63 


123 


72.04 


0.33 


J142154.50+522423.5 


EGS26 


10000. ± 


5464.5 


334.10 ± 129.40 


0.005 


4x10' 


1.75 ± 1.15 


1.50 ±0.50 


63 


4.64 


0.35 


J141717.29+523217.6 


EGS27 


17.4 ± 


59.9 


4.07 ± 


2.72 


0.005 


10 5 


4.00 ± 1.43 


0.00 ±0.13 


400 


228.35 


0.46 



NOTE. — "The redshifts are GALEX-derived spectroscopic redshifts from Deharveng et al. (2008). Coordinates and best-fit modeling results for our 30 
LAEs. The reported 1 a errors are discussed in §2.1. In a few cases the error is formally zero, as these grids are not quite continuous, thus all the middle 68% of 
the probability distribution all fell at the exact same age or dust value. 



we compute a grid of models, varying the metallicity, star for- 
mation history (SFH), stellar population age, and dust extinc- 
tion (using the dust law of Calzetti et al. 1994). The metal- 
licities were allowed to be the six discrete values available in 
BC03, from 0.005 - 2.5Z Q , while the SFHs were exponen- 
tially decaying, with five values of the characteristic decay 
time tsfh ranging from 10 5 - 4 x 10 9 yr. The age grid was 
composed of 62 out of the 244 possible ages in BC03, ranging 
from 1 Myr, up to the age of the Universe at z ~ 0.3 (~ 10 
Gyr). Lya and Ha emission are included using the number 
of ionizing photons (output by BC03) and case B recombi- 
nation. Various levels of interstellar medium dumpiness (i.e. 
Neufeld 1991; Hansen & Oh 2006; Finkelstein et al. 2008, 
2009) are allowed by changing the effective dust optical depth 
for Lya relative to the continuum. We do this by multiplying 
the continuum flux by e" T and the Lya flux by e" f/r , where 
r is the optical depth due to dust, and q specifies the level 
of ISM dumpiness (where q = simulates extremely clumpy, 
and q = 10 simulates homogeneous). This effectively decou- 
ples the Lya line strength from the SED fitting, such that the 
line does not force the best-fit population to be younger than 
the continuum would indicate. 

The model spectra are then redshifted to the GALEX spec- 
troscopic redshift of a given object (in redshift bins of 0.02), 
and attenuated by intergalactic medium (IGM) absorption via 
the Madau (1995) prescription. The best- fit model is found 



via x 2 minimization among the five free parameters, using the 
GALEX Lya line fluxes, FUV and NUV fluxes for all objects, 
along with the optical broadband data discussed above. The 
model fitting is independent of the mass, as we fit flux ratios 
(f A /f,) to find the best-fit model, and then derive the mass by 
the weighted mean of the observed optical fluxes to the best- 
fit model fluxes. See Finkelstein et al. (2008, 2009) for more 
extensive details on the model fitting process. 

In order to estimate uncertainties in our results, we followed 
the method prescribed by Kauffmann et al. (2003), which uti- 
lizes Bayesian likelihoods to estimate model-fitting parame- 
ter uncertainties. The most probable value of a parameter is 

thus the peak of the probability distribution, P oc e~ x " . This is 
simple for the one parameter case, but we have five free pa- 
rameters. Thus, to find the collective error in one parameter, 
we marginalize over all other parameters, in effect collapsing 
the errors down into one dimension. This equates to assuming 
a uniform prior probability distribution in these other param- 
eters. There is no requirement on the shape of this proba- 
bility distribution, thus the 1 a error is the range of parame- 
ters once the 16% tails on each end are removed. Figure 3a 
shows the age probability distribution of one of our objects. 
This method is only useful for those parameters which have a 
quasi-uniform distribution, thus we are restricted to obtaining 
parameter errors to the stellar population age, dust extinction 
and q. We also report errors on our derived masses, how- 
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FIG. 2. — Best-fit single population models to the z ~ 0.3 LAEs, with their observed SEDs over-plotted. The stellar population modeling procedure is discussed 
in §2.1. The open circles represent the bandpass-averaged fluxes of the models, while the asterisks represent the observed data points. We fit the models to the 
GALEX derived Lya line, NUV and FUV fluxes, as well as optical fluxes from the MUSYC and CFHTLS surveys for the ECDF-S and EGS objects, respectively. 
The emission lines shown in the models are Lya and Ha. 
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ever given that our best-fit masses came from fitting the col- 
ors of the objects, we have no mass \ 2 distribution. Rather, 
in order to ensure that the uncertainty in all of the parame- 
ters comes through to the uncertainty in the mass, we ran 10 3 
Monte Carlo simulations, varying the flux by an amount pro- 
portional to the flux error (inclusive of the zeropoint error; see 
§3). The quoted mass error is the range of the middle 68% of 
the 10 3 best-fit masses. 

Prior to model fitting, it is important to ensure that the fluxes 
being used are consistent. Typically, one would prefer to use 
the total flux, which can be approximated by the SExtractor 
MAG_AUTO extraction. Both the CFHTLS and MUSYC 
fluxes used in the paper come from the MAG_AUTO mea- 
surement. While the GALEX FUV and NUV fluxes have 
a photometric accuracy of 0.05 and 0.03 mag, respectively 
(Morrissey et al. 2007), these measurements can be highly un- 
certain given the broad FWHM of the GALEX beam. We thus 
assume a standard error of 0.2 mag for both the GALEX FUV 
and NUV fluxes, which should result in their being assigned a 
lesser weight during the fitting process, as befits their more 
uncertain nature. In addition, while we expect the Madau 
IGM prescription to be valid at these redshifts, if the IGM 
is different at z ~ 0.3, this would affect primarily only the 
FUV flux, showing a further need to reduce its weight in the 
fit. As we mention above, we also use the GALEX measured 
Lya line flux measurement in our model fitting, which we 
believe is well calibrated due to the high throughput of the 
GALEX prism. D08 measure a 1 a precision on the line-flux 
measurement of x 10~ 16 erg s" 1 cm" 2 , which we use during 
our fitting process. 

3. RESULTS 

Table 1 tabulates and Figure 2 displays the best-fit models 
for each of our objects. We gave our objects new identifiers, 
with CDFS1-3 representing the three LAEs in the ECDF-S, 
and EGS1-27 the 27 LAEs in the EGS. While we also fit a 
model allowing two separate bursts of star formation (one at 
10 Gyr ago, and one at any time; see Finkelstein et al. 2009 
for details), these fits were ruled out by a large margin in most 
cases, only providing a marginally better fit to two objects 
(CDFS3 and EGS1) 2 . Since we are mainly interested in the 
average properties of these LAEs, we thus move forward us- 
ing only the single population models. 

During our initial model fitting, the lowest reduced \ 2 (x?) 
was still much greater than one, even though many of the best- 
fit models appeared to match the observations. We found that 
our flux errors were likely underestimated, as these objects 
are so bright, that their flux errors are rather small. To cor- 
rect for this, we added a systematic error equal to the uncer- 
tainty in the photometric zeropoint for each band into the x 2 
denominator, which allowed for the objects with the best fits 
to have x 2 near to unity (e.g., Papovich et al. 2001). These 
errors were 1.2, 0.9, 0.4, 0.5, 0.6 and 1.4% for the MUSYC 
U, B, V, R, I and z' bands, respectively, and 2% for all of 
the CFHTLS bands. We attempted to rule out the possibility 
that any of these objects were AGN by examining X-ray data 
where available (AEGIS Chandra data) to see if any of our ob- 
jects were detected. We found one object which had a nearby 
X-ray counterpart; EGS9, with an X-ray source within 0.2". 

2 In the limit of no mass in the old population, these models are not equal 
to the single-population models, as these models are fixed to have a "burst" 
SFH, which we model with t^fh = 10 5 Y 1 ' 



However, the full-band X-ray luminosity of this source is hps 
= 8.4 x 10 41 erg s -1 , which could indicate a low-luminosity 
AGN, but it is also consistent with X-ray levels from star- 
forming galaxies (e.g., Szokoly et al. 2004), thus we leave this 
object in our sample. It is worth noting that a few other galax- 
ies show possible signs of AGN. EGS4 exhibits a very steep, 
red spectral slope, reminiscent of a power-law. By itself, this 
could just be indicative of heavy reddening. However, the 
inability of the model to fit the GALEX NUV flux could in- 
dicate a possible emission line at that wavelength, possibly 
Hell A1640 or CIII A1909, both of which could be present 
in an AGN. Thus, EGS4 could be a heavily obscured quasar 
(e.g., Gregg et al. 2002; Kuraszkiewicz et al. 2009), how- 
ever we require follow-up optical spectroscopy to know for 
sure. There are a few other objects where the NUV flux is 
difficult to fit, and nebular continuum emission could also be 
expected to contribute at these wavelengths, although this ef- 
fect is likely only significant in metal-poor galaxies within a 
few Myr of a recent starburst (Schaerer 2002). 

3.1. Stellar Population Age and Mass 

From our model fitting, we find an age range of 4 Myr - 
10 Gyr, with a median age of ~ 0.5 Gyr. Correspondingly, 
we derive stellar masses from 1.7 x 10 7 - 3.3 x 10 11 M Q , 
with the lowest masses corresponding to the youngest ages, as 
we would expect. The best fits of the four youngest objects, 
EGS2, EGS25, EGS27 and CDFS2 are all of poor quality (x 2 r 
> 30). Along with the fact that the next youngest object is 
200 Myr old, this casts doubt on the results from these four 
objects. These four objects also have four of the five highest 
values of x 2 (CDFS1 also has x 2 > 30). Excluding these four 
objects from the rest of our analysis, we find age and mass 
ranges of 203 Myr - 10 Gyr (median ~ 1.8 Gyr) and 1.1 x 
10 9 - 3.3 x 10 11 M (median ~ 7 x 10 9 M Q ), respectively. 
This means that the stars in the oldest galaxies in our sam- 
ple formed as long as 10 Gyr ago, or nearly the age of the 
Universe at these redshifts. 

In Finkelstein et al. (2009) we fit stellar populations to 14 
z ~ 4.5 LAEs, and found an age and mass range of 3 - 500 
Myr and 1.6 x 10 8 - 5.0 x 10 10 M , respectively. The lower 
end of our new results is consistent with our previous work, 
while the upper end represents much more evolved objects, 
with masses up to 6.7 times greater. The median stellar mass 
at z ~ 0.3 is ~ 7.5 times greater than the median mass we 
measured at z ~ 4.5. As for age, the median age at z ~ 4.5 was 
12 Myr (with the oldest stellar population age of 500 Myr), 
while at z ~ 0.3, half of the objects are older than 1 Gyr. 

The error analysis discussed above allows us to examine 
the age likelihood distribution of individual objects, by plot- 
ting a histogram for each object to see how many simulations 
resulted in a given age (Figure 3a). We can thus analyze the 
sample as a whole by making this plot for each object, and 
taking the average, which is what we show in Figure 3b. We 
find that an average LAE at z ~ 0.3 is highly likely to have an 
age > 100 Myr, with with many > 1 Gyr. 

Many other studies at high redshift have derived the physi- 
cal properties of LAEs via similar methods, although the ma- 
jority of these studies use stacking analyses to derive average 
properties of the sample, which can wash out individual re- 
sults. In Gawiser et al. (2006a), from a stack of ~ 18 LAEs 
spectroscopically confirmed to be at z ~ 3.1, they found a 
typical age of 90 Myr, and a typical mass of 5 x 10 8 M . 
Also at z ~ 3.1, Nilsson et al. (2007) find a mass of ~ 
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FIG. 3. — (a) The age probability distribution for one of our objects. The asterisk represents the best-fit model (from minimizing x 2 )> while the error bars 
represent the 1 cr range from the probability curve, (b) The distribution of ages from the probability curves averaged over all 30 objects, showing that LAEs at 
low redshift are likely to be > 200 Myr old. 



10 M from a stacking analysis of 24 LAEs, with an age 
from a few to several hundred Myr (the age was not well con- 
strained). Using Hubble Space Telescope grism spectroscopy 
to discover 10 LAEs at z ~ 5, Pirzkal et al. (2007) found 
ages of only a few million years and very low masses, from 
10 6 - 10 8 M Q . Lastly, Pentericci et al (2009) studied Ly- 
man break galaxies (LBGs) with Lya emission (which we can 
also consider to be LAEs, although the LBG selection criteria 
preferentially selects brighter, and thus more massive galax- 
ies), and found stellar masses from ~5 x 10 8 - 5 ~ 10 
M©. Additionally, they found that while the majority of their 
galaxies were best-fit by young populations, a small fraction 
required a significantly older component. Overall, the litera- 
ture is consistent with what we have found at z ~ 4.5, in that 
high-redshift LAEs appear to be mainly young and low-mass, 
although some small fraction do appear to be more involved. 
The evolution of this fraction with redshift is unknown, and is 
a top priority question for future studies. 

3.2. Metallicity, Dust and ISM dumpiness 

The model grid consisted of six allowed values of metallic- 
ity: 0.005, 0.02, 0.2, 0.4, 1.0 and 2.5 Z©. Only two objects 
(CDFS1 and EGS26) out of 26 was best fit with the lowest 
value of metallicity, with five more being best fit with Z = 
0.02 Z©. The majority of objects (19) were best fit with Z > 
0.2 Z©. In F09, we found that at z ~ 4.5 50% of our objects 
were fit with metallicities of 0.02 Z© or less, thus showing 
possible evidence of metallicity evolution with redshift. This 
evolution does not show up in the dust content, as the mean 
extinction at z ~ 0.3 is A1200 — 1.1 mag, while at z ~ 4.5 it 
was A1200 — 2.5 mag. While the ISM dumpiness (q) spans 
a similar range at low redshift as at z ~ 4.5, at z ~ 4.5, over 
50% of the objects had q < 1. At z ~ 0.3, only nine ob- 
jects have dust enhancement of the Lya EW occurring with 
q < 1, but dust is not extinguishing the EWs of these objects 
much either, as 22/26 objects have q < 1.5. Histograms com- 
paring the age, dust, metallicity and dumpiness differences 
between these two redshifts are given in Figure 4. We com- 
puted the Kolmogorov-Smirnov (KS) statistic, using the IDL 
task kstwo, to explore these histograms further. Using age as 
an example, the KS statistic computes the probability that the 
array of age values at z ~ 4.5 and z ~ 0.3 are drawn from 
the same distribution. This returned a probability of 7.5 x 
10~ 7 , thus there is a very low chance that the ages at the two 



redshifts arise from the same distribution. Similarly, we find 
low KS probabilities of 0.003, 0.102 and 0.007 for A1200, Z 
and q respectively, with the slightly higher probability for Z 
possibly due to the very coarse grid of possible metallicities. 
It is comforting that the shape of Figure 3b, and the first panel 
of Figure 4 are similar, as it shows that our ages accurately 
represent the most probable results. 

While Gawiser et al. (2006a) and others (e.g., Lai et al. 
2008) derive Ay consistent with zero in their studies, many 
recent studies do find dust in their LAEs 3 . Even though their 
objects were extremely young, Pirzkal et al. (2007) found 
dust in up to half of their LAEs, with Ay as high as 0.6, while 
Nilsson et al. (2007) found Ay ~ 0.3. Pentericci et al. (2009) 
also found dust in over half of their objects. These results are 
also consistent with our results at z ~ 4.5, although the un- 
certainties in the dust measurements at high redshift are large 
(given that most of these studies only loosely constrain the 
4000 A break), thus any evolution in the dust content in LAEs 
are uncertain. 

These differences in metallicity and ISM properties at dif- 
ferent redshifts may give us insight into the evolution of these 
objects. We would expect the average stellar population to 
become more metal enriched with time, and we believe that is 
what we are seeing with the metallicity differences between 
z ~ 4.5 and 0.3. The fact that objects at z ~ 4.5 have more 
dust extinction on average is interesting, although it may just 
be a selection effect. The z ~ 0.3 LAEs were selected via 
spectroscopy, thus they were more sensitive to low EW ob- 
jects. However, at z ~ 4.5, we required a significantly strong 
narrowband excess in order to select an object as an LAE, 
which biased us towards higher values of the Lya EWs. As 
we showed in F09, a large portion of our sample had high 
EWs possibly due to clumpy dust enhancement, thus we may 
have been biased towards dusty objects as well. According to 
the best-fit models, only 36% of our z ~ 4.5 LAEs had EW 
< 120 A, whereas as z ~ 0.3, over 55% of our LAEs have 
EWs < 120 A. In fact, for z ~ 0.3, the measured EWs from 
D08 may be more accurate as they come from spectroscopy. 
Using these numbers, all 30 z ~ 0.3 LAEs have EWs < 120 
A, with a median value of 31 A, and a minimum of 11.8 A. 
The change in EW with redshift could also be due to the ap- 
parent evolution of these objects, since EW is a ratio between 
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FIG. 4. — Histograms showing the age, dust, metallicity and ISM dumpiness distributions for this z ~ 0.3 sample (solid lines) and our z ~ 4.5 sample (dashed 
lines). Strong differences are seen in all four physical properties, strongly implying that low-redshift LAEs are different than their high-redshift counterparts. 
We list the results of KS tests done on each of the following distributions, showing the probability that the results at the two redshifts are drawn from the same 
distribution. 



the line flux and continuum flux. If the amount of Lya pho- 
tons being produced is similar, but the continuum is brighter, 
the EW will drop. At high redshift, LAEs are known to be 
relatively small, with r e < 1 kpc (Pirzkal et al. 2007), thus it 
is likely that these objects are producing their Lya photons in 
HII regions in much smaller disks than at low redshift. Even if 
low-redshift LAEs had a similar number of HII regions, they 
would be spread over a much larger disk, diluting the EW as 
the larger disk is producing more continuum light. 

One may then ask, why do we not also see the high EW, 
dust enhanced LAEs at low redshift? There are a few, as dis- 
cussed in the above paragraphs, but the numbers are signifi- 
cantly smaller than those at high redshift. This may indicate 
an evolution in the dust geometry in LAEs, and/or galaxies in 
general. Over time, a galaxy will accumulate more and more 
dust as stars continue to grow old and die. Isolated regions 
may become clumpy due to turbulence in H II regions among 
other things. In the high-redshift LAEs, which may be intrin- 
sically smaller, these local events may be enough to change 
the ISM on a galaxy wide scale, allowing for the large range in 
q we found. However, at z ~ 0.3, these objects may be larger 
in size, thus it might be harder for a whole galaxy to have an 
entirely clumpy or entirely homogeneous ISM, which may be 
responsible for the fact that most z ^ 0.3 LAEs have q near 
1 - 1.5. In addition, external events, such as galaxy mergers, 
can severely disrupt the ISM. Ryan et al. (2008) show that the 
major merger fraction peaks at z ~ 1.3. If the time since their 
last merger is large, this will also allow low-redshift LAEs to 



have more settled ISMs. 

Although not much work has been done on LAEs at low 
redshift, a recent study has been published studying the role of 
dust in local Lya emitting galaxies (Atek et al. 2008). Specif- 
ically, their plots of Lya EW vs. dust extinction allow an 
opportunity to check if clumpy dust is affecting the EW, and 
by how much. In Figure 3 of their paper, they show this plot 
for Haro 1 1 . The emission component of this plot shows sig- 
nificant scatter at low E(B-V), however, it approaches a con- 
stant value of - 20 A from E(B - V) = 0.5 - 1.0. In order 
for the Lya EW to stay constant with increasing dust con- 
tent, q must equal 1 , which is the most likely value we found 
for q at z ~ 0.3. This same effect is seen in their plots of 
IRAS 08339+6517, and also in NGC 6090, albeit less signif- 
icantly. We believe that this confirms our prediction that the 
ISM of lower-redshift LAEs are neither extremely clumpy nor 
entirely homogeneous. 

3.3. Stellar Evolution Uncertainties 

While the BC03 models are very widely used across ex- 
tragalactic astrophysics, recent insights into the contributions 
of post main-sequence evolution to the overall SED of a 
galaxy has called these models into question. The BC03 mod- 
els do include all phases of stellar evolution from the zero- 
age main-sequence until the stellar remnant stage, including 
the thermally-pulsating asymptotic giant branch (TP-AGB) 
phase. However, Maraston et al. (2005, 2006) have recently 
shown that the BC03 models underestimate the level at which 
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TP-AGB stars can dominate the rest-frame near-IR (A > 1 
/im), causing overestimates of the derived ages and masses. 

Even prior to our model fitting we would expect these z ~ 
0.3 LAEs to be older, possibly in the range where TP-AGB 
stars can dominate their NIR (~ 0.2 - 2 Gyr). However, our 
reddest band (z') corresponds to rest-frame ~ 7000 A, thus 
we would not expect a different treatment of TP-AGB stars 
to significantly alter our results. Chariot & Bruzual (2007 in 
prep) have devised new models, including the Maraston et al. 
treatment of TP-AGB stars, which will be released in the near 
future. Future studies at redder wavelengths will need to use 
these, or similar, models to more accurately fit their objects. 

4. CONCLUSIONS 

We have performed stellar population model fitting on a 
sample of 30 GALEX discovered z ~ 0.3 Lya emitting galax- 
ies. We derived the stellar population ages and masses, as well 
as metallicities, dust extinction, and ISM geometries. The age 
and mass results imply that low-redshift LAEs are both older 
and more massive than their high-redshift counterparts, with 
the youngest low-redshift LAE having a best-fit age of ~ 200 
Myr, much older than the average age at high redshift. We 
also found that the metallicity in LAEs is higher at low red- 
shift, showing evolution of the mean metallicity of the stars 



inside these galaxies, as would be expected over a period of 
~ 9 Gyr of cosmic time. 

Perhaps most interesting, we found that the ISM at low red- 
shift appears to be very stable, in a quasi-clumpy state. At 
high redshift, we found that the ISM in LAEs ran the whole 
range from extremely clumpy to very homogeneous. At z ~ 
0.3, we found that while a few galaxies were fairly clumpy, 
most resided in the regime of q ~ 1 , where dust will not af- 
fect the Lya EW. This result is confirmed by the study of 
local LAE analogs by Atek et al (2008). Our results show 
that LAEs at z ~ 0.3 show significant differences from their 
high-redshift cousins, although further study of a larger sam- 
ple will help confirm this. Although low-redshift LAE studies 
are difficult as they need to be done from space, help should 
be on the way with the near-UV sensitive Wide Field Camera 
3 (WFC3) due to be installed on the Hubble Space Telescope. 
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